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Abstract. 

The scope of this talk is to present some preliminary results on an effort, currently in 
progress, to generate an exact solution of Einstein's equation, suitable for describing spacetime 
around a rotating compact object. Specifically, the form of the Ernst potential on the symmetry 
axis and its connection with the multipole moments is discussed thoroughly. The way to 
calculate the multipole moments of spacetime directly from the value of the Ernst potential 
on the symmetry axis is presented. Finally, a mixed ansatz is formed for the Ernst potential 
including parameters additional to the ones dictated by Sibgatullin. Thus, we believe that this 
talk can also serve as a comment on choosing the appropriate ansatz for the Ernst potential. 



1. Introduction 

In order to analytically describe spacetime around a rotating compact object in the framework of 
general relativity, one needs an exact asymptotically flat solution of Einstein-Maxwell equations. 
To simplify the problem it is reasonable to assume that a long-lived, quiescent rotating body 
possesses axial symmetry and reflectional symmetry with respect to the equatorial plane. The 
solution must depend on a number of arbitrary physical parameters such as the mass, the angular 
momentum, the mass-quadrupole moment, the magnetic dipole e.t.c. of the central body. 

Sibgatullin has proposed a way to generate stationary axisymmetric solutions pQ, starting 
from the value of the Ernst potentials [21 13] £ and $ on the symmetry axis. His method has been 
used successfully by many authors (see for example [IJ [3] ) . However the arbitrary parameters 
appearing in the proposed ansatz for £ and $ are physically meaningful only in specific and 
rather special cases (for example the Kerr solution). 

On the other hand, it is well known for quite some time that spacetime geometry around 
an axisymmetric object leads to a specific set of the Geroch-Hansen multipole moments (S1IZI- 
Thus, the scalar moments can be used to characterize the spacetime, in the same way that the 
newtonian gravitational moments characterize the newtonian gravitational field. In this paper 
we will try to combine multipole moments with Sibgatullin's arbitrary parameters in order to 
form an ansatz for the Ernst potential £ on the symmetry axis, suitable to obtain a metric for an 
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axisymmetric compact object through Sibgatullin's method. Due to lack of space and time this 
study will be restricted to vacuum spacetime. Note however that it can be easily generalized to 
include electrovacuum spacetimes. 

The plan of this paper is the following: In Sec. |2]we give an outline of Sibgatullin's method. 
In Sec. 01 we show the way to calculate the moments directly from the value of the Ernst 
potential on the symmetry axis. We also present a way to generate a solution with specific 
multipole moments. In Sec. @] we present a "mixed" ansatz for the Ernst potential including 
both Sibgatullin arbitrary parameter and a number of new parameters related to the moments. 
Finally we comment on the restrictions that come up if we additionally impose reflectional 
symmetry with respect to the equatorial plane. 

2. Sibgatullin's method for generating exact solutions 

It is known for quite some time that in the case of a stationary, axisymetric spacetime the metric 
functions can be fully determined if the so called Ernst potential £ is known Sibgatullin 
has proposed a way of generating solutions of Einstein's field equations, that reduces to 
evaluating the solution of a linear system of algebraic equations. His approach is based on the 
definition of the Ernst potential on the symmetry axis in terms of rational functions, involving 
a series of parameters dj (I] 
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We will merely outline some key points of the method. Sibgatullin has shown that the Ernst 
potential £(p, z) can be evaluated from the potential Q which is given on the symmetry axis, 
as the integral 

l 

f fi(a)e(£)da 

£(p,z)=^- / == r (2) 

were the function /i(cr) satisfies the following linear integral equation and normalization condition 
respectively: 

p/^)(e(0 + g(Q) jj^)d L = 7T 

J (a-T)Vl^ ^VT^ KJ 

were the P denotes the principal value of the integral, ^ = z + zap, a,r £ [—1,1] and the 
functions e(£) and e(£) are analytical continuations on the complex plane for the functions e(z) 
and e*(z). The first step is to express equation Q as an expansion of elementary fractions e(z) = 
1 + Sj=i an d the sought for solution of the integral equation as p(a) = A)+Z)l=i J~§^ were 
the dj and are the roots of the polynomials z n + Mz"" 1 + Y^j=i o,jZ n ~^ = and e(z)e*(z) = 
respectivly . Substituting the series expansions of e(z) and ju(cr) in the integral equation and the 
normalization condition we obtain a closed system of 2n + 1 linear algebraic equations for the 
2n + 1 parameters A}.. These coefficients can be evaluated and they completely determine the 
Ernst potential that will generate the solution. In particular, it has been shown by Ruiz, Manko 
and Martin [HJ that the metric functions can be evaluated with the use of some determinants 
which are functions of the parameters e,- , dj , . 

3. The multipole moments of the solutions 

In |Hj an algorithm is presented for the evaluation of the multipole moments of stationary 
axisymmetric spacetimes from the series expasion of the function £ on the symmetry axis. The 



real and the imaginary part of £ are related to the so-called Hansen potentials respectively jjj. 
Since the multipole moments are evaluated at infinity, we will use the prefered coordinates 

(4) 
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where p and z are Weyl coordinates. As mentioned in (Hj £ can expanded around z = in the 
following way. 

oo 

i(p = 0) = J2mnZ n (5) 
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The multipole moments can then be evaluated from the m n 's using the relations presented in 
p. 2255 of 0. 

If we express £ in terms of the Ernst potential £ then 

,-, „> i-£(p = o) 
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or in the prefered coordinates 



It is now obvious that the moments can be evaluated from <?(/? = 0). Using (jSJ) and Q we get 

$>^ = . (8) 
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Equation can be used to relate nii with and vise versa. Applying Taylor expansion at 
z = we derive the following recursive formulae: 



m = M, m n = - a l m n~l or a = 1, a n = - ^ T7 a n-«- (9) 



To derive the second set of expressions , of course, one has to regard that the number of m n 's is 
finite and the number of a ra 's is infinite. Before we make any comments we present two examples. 
If we set 

a n = 0, n = 1,2,3,... (10) 

then it follows that 

m = M, m n = 0, n= 1,2,3,... (11) 
which will lead to the Schwarzchild solution, while if we set 

a\ = — ia, a n = 0, n = 2,3,4,... (12) 

then it follows that 

M 2n = (-l) n Ma 2n , J 2n+1 = (-l) n Ma 2n+ \ n = 0, 1, 2, 3, . . . (13) 
which are the mass and the mass current moments of the Kerr solution respectivelly. 



We have related £(p = 0) and thus Sibgatullin's arbitrary parameters with the moments. 
This is another way to demonstrate that all information regarding spacetime are included in the 
value of the Ernst potential on the symmetry axis. Using eqs. ® one can compute the moments 
without much effort even before generating the solution. One can also get an insight for the 
number and the nature of the a, 's needed to generate a solution with a certain set of multipoles. 
However, it is important to say that a physical object is in principle expected to have an infinite 
number of moments. It is easy to see that a finite nunber of Oj's imposes a relation between 
the different order multipoles. The more aj's one uses the more flexible this relation becomes. 
However every new aj included leads to an increace of the order of the polynomials of the Ernst 
potential, and consequently creates serious difficulties in generating the corresponding metric. 

Before we close this section it is important to notice that Sibgatullin's method can be also 
used to generate an exact solution with a desired set of multipole moments or equivelantly a 
desired set of mi's. If we solve J7J) for £(p = 0) we get 
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We can use (|T1|) instead of Q and follow the procedure described in section El to obtain the 
solution. We could turn the infinite sum of eq. ® into a finite one, by setting all the moments 
of order > n equal to zero. This, even though not exact, can most of the times be considered as 
a good approximation. 



4. A mixed ansatz for the Ernst Potential 

Sibgatullin's ansatz for £(p = 0) is rather general, but not the easiest one can use to generate a 
solution as we will show. Consider the following ansatz for £ (p = 0): 

n n 
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It is easy to see that if aj = for all j's then bi = rrii and if all bi = the ansatz is reduced to 
Sibgatullin's one. It is obvious that it is not possible to define a new finite set of parameters a'j 
in order to bring the ansatz in the form of eq. (^Q) . This simply depicts the close relation of the 
frj's with the m^s. Notice that one can use the above ansatz to generate an exact solution with 
Sibgatullin's method. No modifications in the generation algorithm are necessary. The merit of 
using such a mixed ansatz with I a's and n 6's lies in the fact that the degree of its polynomials 
will be equal either to / or to n + 1 (depending on which one is the largest), whereas in the 
original ansatz it is equal to k = I + n in order for it to have the same number of parameters. 
This of course can lead to serious simplification of the computations needed to generate the 
solution. For example a 6 parameter solution would require third degree polynomial instead of 
fifth. 

Before we proceed lets take a moment to focus on the special case in which spacetime has 
an additional symmetry: reflection symmetry with respect to the equatorial plane. This is a 
realistic assumption for most astrophysical objects. In this case rrii should be real for even and 
imaginary for odd i's CHI EH- By evaluating the mj's in terms of the a^'s and b^s one can 
show that, if there is a finite number of a^'s and ftj's then they should also be real for even i's 
and imaginary for odd i's. 



A simple case for an ansatz is the following: 



_ , _ z 2 - Mz - iaz + bc- icM 
MP-U)- z 2 +Mz _ iaz + bc + icM > W 

where all parameters are real. The first five multipole moments of the corresponding spacetime 
are the following: 

M = M, Mi = 0, M 2 = -(a 2 + ac + bc)M, M 3 = 0, 

M 4 = (a 4 + a 3 c + 3a 2 bc + 2abc 2 + 6 2 c 2 ) M ~-=[m ((a + c) 2 M 2 - (a 2 + ac + 6c) M 2 ) 
J = 0, Ji = (a + c)M, J 2 = 0, 

J 3 = -a 3 M - a 2 cM - 2abcM - bc 2 M, J 4 = 0. (17) 

If we set c = we get 

M 2n = (-l) n Ma 2n , J 2n+1 = (-l) n Ma 2n+ \ n = 0, 1,2,3, ... . (18) 

On the other hand if we set b = a + c 

M 2n = (-l) n M(a + c) 2n , J 2n+1 = (-l) n M(a + c) 2n+ \ n = 0, 1, 2, 3, . . . . (19) 

Therefore, there are two discrete ways for the metric to be reduced to the Kerr metric, having 
different angular momentum. The first one is by reducing the ansatz to the known ansatz for 
the Kerr metric by setting c = as we showed. The second one which does also reduce to a 
Kerr metric with angular momentum c if a = 0, is by setting b = a + c. From eq. ((17)) one 
can see that by using freely the parameters a, 6, c he can give the desired value to the angular 
momentum and the mass quadropole and still have two free parameters including the mass. 
Thus we think that starting from this ansatz it will be possible to generate a solution suitable 
to describe spacetime around compact rotating objects sufficiently well. 

Finally, it is worth mentioning that this whole approach can be generalized to include 
electrovacuum spacetimes (based on |31 fTTl IT!?] ). 
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